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Super-oscillating signals are band limited signals that oscillate in some region faster their largest 
Fourier component. While such signals have many scientific and technological applications, their actual 



use is hampered by the fact that an overwhelming proportion of the energy goes into that part of the 
signal, which is not super-oscillating. In the present article we consider the problem of optimization of 
such signals. The optimization that we describe here is that of the super-oscillation yield, the ratio of 
the energy in the super-oscillations to the total energy of the signal, given the range and frequency of 
the super-oscillations. The constrained optimization leads to a generalized eigenvalue problem, which is 
solved numerically. It is noteworthy that it is possible, to still increase the super-oscillation yield at the 
cost of slightly deforming the oscillatory part of the signal, while keeping the average frequency. We 
^sO | show, how this can be done gradually, which enables a trade-off between the distortion and the yield. 
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Super oscillatory functions provide a stunning refutation of a very widely accepted lore, that band 
limited functions cannot oscillate with a frequency larger than its maximal Fourier component. A number 
of examples have been given in the past for such functions with very interesting applications to Quantum 
Mechanics CD-El, signal processing flTI- lfTTI and to optics, where super-oscillations are intimately 
related to super-resolution fl2l - fl7l . Interestingly, it was discovered that in random functions, defined 
as superpositions of plane waves with random complex amplitudes and directions, considerable regions 
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are naturally super-oscillatory lTT8l . flT9l . Various mathematical aspects of the phenomenon have been 
discussed also more recently in EOl . |[2TI . 

In a very important sense, though, the idea that a band limited function cannot oscillate faster than 
its largest Fourier component is not entirely false. It is well known that the super-oscillations exist in 
limited intervals of time (or regions of space, depending on the actual problem) and that the amplitude 
of the super-oscillations in those regions is extremely small compared to typical values of the amplitude 
in non-oscillating regions ||2), ifTUl . ifTTTl . It is generically so small, that any hope of practical application 
of super-oscillating functions depends on tailoring the functions carefully to reduce that effect as much 
as possible. Two different approaches have been offered over the years to the problem of optimization 
of super-oscillation Q, IfTTTl . which will discussed in relation to ours as we proceed. 

In contrast to the two last references we prefer to work with periodic functions, which will thus be 
described by a finite number of Fourier coefficients. This choice, we believe has a number of advantages. 
First, it is clear, that it is more difficult to achieve super-oscillations, with a finite number of degrees 
of freedom than with an infinite number of degrees of freedom encoded in the Fourier transform of a 
band limited non-periodic function. Thus, achieving super-oscillations with a finite number of degrees 
of freedom is more challenging, yet enables a clear view of questions concerning the total number of 
oscillation versus the number of degrees of freedom. Second, we obtain an easy and practical way of 
constructing optimal super-oscillations. 

Consider the function 

/(t) = -^+E^coB(mt) (D 

V27T ~! V* 

Choose an interval [0, a] with a < 7r and impose on the function M constraints in the interval, f(tj) = fj,j 
for < tj < a and j = 0, . . . , M — 1. The constraints result in a set of M linear equations in N + 1 
unknowns of the form, 

TV 

^ ] Cj m A m = Cj • A = //j , (2) 

m=0 



/IT 



where 

( cos(mtj), for m/0 
l/v2, for m = 

Generically this set of equations has no solution for M > N + 1, has one solution for M = N + 1 
and a whole space of solutions for M < N + 1. In particular, we can choose 
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Provided M < N + 1, this choice constrains the function to oscillate within the interval [—a, a] between 
the values ±1 with a frequency 

M= '' Jtf - 1 '. (4) 
a 

It is thus clear that the frequency of oscillation within the interval [—a, a] can be increased indefinitely 
just by decreasing its size. Therefore, although to have a solution at all, we need that M < N + 1, the 
ratio between uj and N, the largest frequency appearing in the Fourier series, can be made as large as 
we want. Thus it is not a problem at all to obtain super-oscillations. This comes, at a cost, of course. 
First, we can obtain super-oscillations with a prescribed frequency ui only within an interval [—a, a] with 
a < tlK an d as stated before and will be demonstrated in the following (see Fig. Q] below) the amplitude 
in that region is relatively extremely small. 

Next, we would like to optimize our super-oscillating function for fixed a and M < N + 1 but we 
have to decide first in what sense do we want to optimize it. Ferreira and Kempf [11], consider the 

oo 

energy of the signal, E = J f 2 (t)dt, use the fact that / is band limited and minimize the energy under 

— oo 

the interpolation constraints (equation ©). We believe that for many applications, the right quantity to 
maximize under the constraints is the super-oscillation yield, 

/ f 2 (t)dt 

Y(M,a) = -£ , (5) 

I P(t)dt 



rather than the total energy. (Note that as will become evident in the following the yield is not just a 
function of uj but of M and a separately. For the discrete case described in (Q]), we take instead of the 
energy, which is infinite, the energy per period. Thus the super-oscillation yield that we maximize under 
the constraints is 

a N 

I f{t)dt £ A mn A m A k 
Y(N, M, a) = ~4 = ^ = L, (6) 

— 7T m=0 



where the entries of the matrix A are given by 

' 2"icos(ria) sin(ma)-n cos(ma) sin(na) 



1 / sin(2na)\ 
' a H ) , m = n / 



m = n = 
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— sin(na), m = 0, nf 

nn 

a 

; 

7T 



(V) 



October 9, 2012 



DRAFT 



4 



The set of M vectors {Cj} defined in Q spans an M dimensional space. In this space we introduce 
an orthonormal basis {ejv-M+i, • • • > &n}- An orthonormal basis for the full N dimensional vector space 
is then constructed by adding the set {eo, • • • , such that e»j • ej = S^j for all i, j = 0, . . . , iV (up 

to the orthogonality requirement, the basis vectors can be chosen randomly, of course). We obtain thus 
the rotated degrees of freedom, 

Bi = ei-A, (8) 

with the obvious advantage that the last M degrees of freedom in B are constrained independently of 
each other and are equal to linear combinations of the /Xj-'s. Thus we denote 

Bi = [ii for i = N-M + l,...,N. (9) 

The numerator / in Q can be written now in terms of the rotated degrees of freedom B as / = 

N N 

E A mn A m A n = J2 &mn R >B m B n , where A (R) = RAR 1 , R being the rotation that takes A 

m.n=0 m,n=0 

into B. Let us describe the matrix A^ R ) by the following block form 



^ A(N+l-M)x(N+l-M) ^(N+l-M)xM 
y Fmx(N+X-M) A-MxM 



(10) 



where V is the transpose of T. The super-oscillation yield expressed in terms of the unconstrained B's 
is 

N-M N-M N N 

m n=0 m = Tl= m,n= 

v _ N+l-M N+l-M n1 x 

tjt = m=N + l — M 

Differentiating the yield with respect to B m and equating to zero yields 

B = -(A-YI)- 1 r/2, (12) 

where I is the unit matrix and jl is the vector of the /L,'s. It is thus clear, that the components of the 
vector B depend on / and on D only through the ratio Y = j^. Those components are, by Cramers 
rule, a ratio of two determinants. The determinant in the denominator is clearly a polynomial of degree 
N + 1 — M in Y. For each entry of B the determinant in the numerator is that of the matrix, obtained 
from (A — YT) by replacing one of the columns by the vector (Tjl). Therefore, B m , the m's component 
of the vector B, for which the yield is extremal, is given by the explicit expression 

pN-M (y\ 
R — m ll I 

D m - p N+1 - M (Y) ' K } 
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Fig. 1. The super-oscillating signal for N = 10, M = 5 and a — 1 (a) on a linear scale in (— 7T, 7r), (b) on a linear scale 
focusing on the segment ( — 1, 1), and (c) on a semi-logarithmic scale. 



where Pj^~ M (Y) is a polynomial of degree N — M in Y and p N + 1 ~ M (y) is a polynomial of degree 
N + 1 — M. Clearly, an equation for the super-oscillation yield can be obtained by plugging the right 
hand side of (fl2l ) directly into the right hand side of (fTOb . We prefer, though, a different rout that yields 
a more simplified form for the equation. The expression C = D 2 B m 9 q^ D ^ is identically zero by 

m=0 

the extremum condition. This results in a simplified equation for Y, 

N-M M / M \ 

E E r™/^~ M cn = E K] yp n ^- m (y) 

m=0 n=l \n=l / 

M \ 

E A mn AmAn P N +l~ M (Y) , (14) 
im,n=l / 

which is a polynomial equation of degree N + 2 — M. We view each of the solution as a generalized 
eigenvalue, as to each of the roots corresponds a N + 1 — M dimensional vector. (In contrast to the 
traditional eigenvalue problem, the eigenvectors are determined by the inhomogeneous linear equation 
(fTTT) and the number of generalized eigenvalues is N+2 — M, and those determine N + 2 — M generalized 
eigenvectors. This means, of course, that the set of generalized eigenvectors is linearly dependent.) We 
are naturally interested in the largest generalized eigenvalue that corresponds to the maximal super- 
oscillating yield. The generalized eigenvector corresponding to that solution has in the interval [—a, a] 
the exact super-oscillation frequency imposed by the constraints. In Fig. Q] we present the super-oscillating 
signal for iV = 10, M = 5 and a = 1. 

A natural question to ask is: can we learn anything from generalized eigenvectors corresponding to 
lower generalized eigenvalues? In Fig. [2] we present the super-oscillating portion of the signal correspond- 
ing to various generalized eigenvalues for N = 10, M = 6 and a = 1. It is obvious that as we go to 
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Fig. 2. The super-oscillating portion of signals corresponding to the various generalized eigenvalues, with Ae (figure (a)) being 
the maximal eigenvalue, which is the maximal yield and Ai (figure (f)) the smallest. The red circles represent the constrained 
points. The eigenvalues are (a) A 6 = 0.00233, (b) A 5 = 4.847 x 10~ 7 , (c) A 4 = 1.379 x 10" 10 , (d) A 3 = 2.559 x 10~ 14 , (e) 
A 2 = 2.176 x 10~ 18 , (f) Ai = 1.189 x 10~ 23 . As can be seen, every time we switch to a smaller eigenvalue another oscillation 
inside the super-oscillating region appears. 

lower generalized eigenvalues the number of oscillations in the super-oscillating interval keeps increasing, 
adding more oscillations to those imposed by the constraints. In fact the number of oscillations in the 
interval [—a, a] grows exactly by one when we go from a generalized eigenvalue Aj to the one immediately 
below it Aj_i < Aj. The number of oscillations inside [—a, a] corresponding to the generalized eigenvalue 
Ai is thus exactly N + 1 - % (1 < i < N + 2 - M). 

How is \i(N, M, a) related to the maximal generalized eigenvalue, Y(N, M, a), for the case with the 
same N and a but with a constrained number of oscillations equal to the actual number of oscillation 
corresponding to Aj, N + 1 — i7 The answer seems intuitively clear. Obtaining the same number of 
oscillations, while imposing less constraints, is expected to give a higher yield. Namely, we expect that 
for every M' > M and 1 < i < N + 2 - M' it will hold that A,(iV, M, a) > Xi(N, M', a). In particular 
since Y(N, M, a) = Xn+2-m(N, M, a) we also have Aj(iV, M, a) > Y(N, N + 2 - i, a). This intuitive 
feeling, is obviously exact, in the case where the set of M constraints, S(M) obeys S(M) C S(N+2 — i). 
To complete the picture for cases where S(M) <£_ S(N + 2 — i), we give more details in Figs. [3] and 
@1 In Fig. [3] we give the different eigenvalues for fixed N and M as a function of a. In Fig. [4] we give 
the eigenvalue for fixed a and N as a function of i for various values of M . This will give not just an 
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Fig. 3. The dependence of the eigenvalue Xi on a for a < 0.5, fixed N — 10 and M = 5. 
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Fig. 4. The dependence of the eigenvalue Xi on the index i for fixed N — 10, a = 1/64 and various M's. 



inequality but a full quantitative picture for small (not necessarily very small) a. 

The improvement in the yield for a given number of oscillations in the super-oscillation region, obtained 
by decreasing i comes at a price, though. The seemingly periodic structure near the middle of the super- 
oscillation interval, corresponding to the highest eigenvalue, generically deteriorates as the number of 
super-oscillations increases. This suggests a possible trade-off between the quality of super-oscillation 
and the yield. For that trade-off, we have at our disposal a whole spectrum of signals, with the same 
number of super-oscillations, corresponding to the same number of low Fourier components. On one 
side of the spectrum we find the signal where all the oscillations are constrained, such that the quality 
of the super-oscillation is good but the yield is relatively low. On the other side of the spectrum we 
find the function with M = 1. For that function we only fix the value of / at the origin to be 1. This 
imposes no oscillation at all on the function in the super-oscillating region. Since, optimal f's, which are 
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necessarily symmetric, are not expected to vanish at the origin, this implies that we are not constraining 
the function at all and the maximization of the yield is equivalent to maximizing it under the requirement 
that the total energy is normalized. This will yield a set of N + 1 ordinary eigenvalues and eigenvectors. 
In the corresponding continuum problem the number of eigenvalues is infinite and the eigenvectors are 
the prolate spheroidal wave functions of Slepian and Pollak |7]. (Ref. |9] mentions a discrete case but it 
is a totally different discreteness than that we study, i.e., in equation (Q]).) It is interesting to note that the 
signals obtained by the discrete Ferreira-Kempf procedure ifTTi do not belong to the family of functions 
described above. Those signals are obtained by minimizing our denominator (defined in ©), 

N N—M N 

D=J2^m= E B l+ E /&> (is) 

m=0 m=0 m=N+\—M 

and under the same oscillation constraints we use (i.e., equation (0). It is clear thus that those eigenvectors 
are obtained by simply setting 

B m = for m = 0, ...,N-M. (16) 

To conclude, we have shown in this paper how to obtain optimal super-oscillating signals that allow a 
gradual trade-off between super-oscillation yield and quality of the signal. Since our optimization process 
is based on a specific way of constraining the signal to produce super-oscillations, given by equation ([3]>, 
improvements of the yield and/or the super-oscillation quality may be expected and will be discussed in 
future publications. 
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